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Abstract — Fast resonance frequency modulation of a supercon- 
ducting stripline resonator is investigated. The experiments are 
performed using a novel device which integrates a hot electron 
detector (HED) into a superconducting stripline ring resonator. 
Frequency modulation is demonstrated by both applying dc 
current or voltage to the HED, and by applying optical illumi- 
nation, with modulation frequencies of up to 4.2GHz. Potential 
applications for such a device are in telecommunication, quantum 
cryptography and biofluorescence. 

Index Terms — Parametric excitation, superconducting ring 
resonator, hot-electron detector, optoelectronics, NbN. 

I. Introduction 

RESONANCE parametric amplifiers are characterized by 
very low noise, high gain, and phase sensitive ampli- 
fication. Parametric resonance in superconducting resonators 
[1] may allow some intriguing applications such as quantum 
squeezing [2], quantum non-demolition measurements [3], 
photon creation by the so-called dynamical Casimir effect [4], 
and more. 

Parametric excitation occurs when the resonance frequency 
of an oscillator varies in time. The first parametric res- 
onance occurs when the excitation is performed periodi- 
cally at twice the resonance frequency fo, namely /(<) = 
fo [1 + £ cos (47r/oi)] [5]. The system's response to such an 
excitation depends on the dimensionless parameter £Q , where 
Q is the quality factor of the resonator. When £Q < 1 the 
system is said to be in the subthreshold regime, while above 
threshold, when £Q > 1, the system breaks into oscillation. 
Achieving parametric gain where £Q > 1 requires that the 
shift in the resonance frequency exceeds the width of its peak 
[6]. 

The frequency modulation mechanism, we employ here, is 
based on changing the boundary conditions of a superconduct- 
ing resonator. This is done by switching a small section of the 
resonator to a normal state by using optical illumination. The 
■ switching time in superconductors is usually limited by the 
relaxation process of high-energy quasi-particles, also called 
'hot-electrons', giving their energy to the lattice, and recom- 
bining to form Cooper pairs. Recent experiments with pho- 
todetectors, based on a thin layer of superconducting Niobium- 
Nitride (NbN), have demonstrated an intrinsic switching time 
on the order of 30ps and a counting rate exceeding 2GHz 
(see [7] and references therein). Resonance frequency shift 
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by optical radiation [8], [9], [10], or high-energy particles 
[11], [12] (for which the required condition £Q = 1 has been 
achieved) was demonstrated, though no periodic modulation 
was reported. Resonance frequency tuning [13] and switching 
[14] as well as optical and microwave signal mixing [15], [16] 
were demonstrated in normal-conducting GaAs microstrip ring 
resonators. 

In this paper we show experimentally that resonance fre- 
quency modulation, at twice the resonance frequency, is within 
reach using superconducting microwave resonators. Further- 
more, the parametric gain threshold conditions, namely £,Q > 
1, is demonstrated in a continuous wave (CW) measurement. 
The experiments are performed using a novel device, that 
integrates a HED into a superconducting ring resonator. The 
HED is used as an optically tuned, lumped element, that 
changes the boundary conditions of the resonator [17], and 
thus manipulates its resonance frequencies. 

In the following section we describe the circuit design and 
fabrication process. The results section starts with the HED 
response to applied dc voltage and current, and the resulting 
effect on the resonance frequencies. These results are followed 
by a comparison with a theoretical model. Afterwards, the 
effect of CW and modulated infrared (IR) light on the reso- 
nance frequency is described, and fast optical modulation of 
the resonance frequency is demonstrated. 

II. Circuit Design and Fabrication 
A. Circuit Design 

The circuit layout is illustrated in Fig. [JJa). The device is 
made of 8nm thick NbN stripline, fabricated on a sapphire 
wafer, with dimensions of 34 x 30 x 1mm 3 . The design 
integrates three components. The first is a superconducting 
ring resonator and its feedline. Ring configuration is a sym- 
metric and compact geometry, which is generally suitable for 
applications, which require resonance tuning [13]. The first 
few resonance frequencies are designed for the S&C bands 
(2 — 8GHz). The resonator is weakly coupled to its feedline, 
where the coupling gap is 0.4mm. The stripline width is set 
to 347/zm, to obtain a characteristic impedance of Zq — 50 ft. 

The second component is a HED, which is monolithic ally 
integrated into the ring structure. Its angular location, relative 
to the feedline coupling location, maximizes the RF current 
amplitude flowing through it, and thus maximizes its coupling 
to the resonator. The HED, shown in Fig. [Tib), has a 4 x 4/im 2 
meander structure, consists of nine NbN superconducting 
strips. Each strip has a characteristic area of 0.15 x 4/.tm 2 and 
the strips are separated one from another by approximately 
0.25/um [18]. 
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Fig. 1. (a) Device layout, (b) Optical microscope image of the HED. 



The HED operating point is maintained by applying dc bias. 
The dc bias lines, forming the third component, are designed 
as two superconducting on-chip low-pass filters (LPF) with 
a cut-off frequency of 1.2GHz. As this frequency is lower 
than the fundamental resonance frequency of the resonator, the 
intrinsic fields of the resonator are not appreciably perturbed. 
A cut of 20/im is made in the perimeter of the resonator, to 
force the dc bias current flow through the HED. 

The device, which is top covered by a bare sapphire 
substrate, is housed in a gold plated Faraday package made of 
Oxygen Free High Conductivity (OFHC) Copper. Supercon- 
ducting Niobium ground planes are dc-magnetron sputtered on 
the inner covers of the package. RF power is fed using a SMA 
launcher, coupled to the feedline. A dc bias is fed through two 
7r-LPFs, screwed to the package, having a cut-off frequency of 
1MHz. IR laser light is guided to the device by a fiber optic 
cable. A through hole of 1mm in diameter, is drilled in the 
Faraday package, and a fiber optic connector affixes the tip of 
the fiber cable at approximately 9.55mm above the HED. 

B. Fabrication Process 

The fabrication process starts with a thorough pre-cleaning 
of the sapphire wafer in solvents. We have experienced that 
the commonly employed process of piranha followed by 
RCA cleaning substantially reduces the NbN adhesion to the 
Sapphire wafer. In the next step, 200nm thick gold pads are 
thermally evaporated through a mechanical mask to form the 
dc contact pads. The mask partially shadows the evaporation 
and thus the pads' perimeters are smoothed. Epitaxy, 8nm 
thick, NbN film is then deposited at 700° C using a dc- 
magnetron sputtering system [19]. Sputtering parameters are 
summarized in table HI and the process itself is further detailed 
in [20]. Next, an A1N layer of 7nm thickness is in-situ 



TABLE I 
Sputtering Parameters 



Process parameter 


NbN 


A1N 


Partial flow ratios (Ar,N2) 


(87.5%,12.5%) 


(0%,100%) 


Base temperature 


700 °C 


60 °C 


Base pressure 


3.6 • 10~ 7 ton 


1.8 • 10" 7 torr 


Work pressure 


6.8 ■ 10~ 3 torr 


2.9 ■ 10" 3 torr 


Discharge current 


360mA 


360mA 


Discharge voltage 


311V 


292V 


Deposition rate 


3.6A/sec 


l.lA/sec 


Thickness (t) 


8nm 


7nm 


Target-substrate distance 


95mm 


230mm 



sputtered in N2 atmosphere at a temperature smaller than 
100°C. This layer protects the vulnerable NbN layer during the 
following fabrication processes and restrain degradation [21]. 
It has also a functional role, as at cryogenic temperatures, 
it serves as a thermal conducting layer, which enhances the 
cooling of the NbN layer. In the next step, the HED meander 
stripline is patterned using electron beam lithography (EBL). 
The deposition of a 80nm thick PMMA 950K layer is followed 
by EBL with the following parameters: 40kV, 15pA, and 
2.1nC/cm, corresponding to acceleration voltage, emission 
current, and line dose respectively. Afterwards, the A1N layer 
is directly etched through the PMMA mask using ion milling. 
The remaining A1N layer serves as a mask for the sequential 
etching of the NbN layer, using low power reactive ion etching 
(RIE) in SFg environment [22]. The remaining PMMA is 
removed by NMP The last fabrication step is the patterning of 
the resonator and the LPFs features. This is achieved by using 
standard photolithography process. The photoresist develop- 
ment process (employing AZ-326 photoresist developer), also 
wet etches the A1N layer, while the remaining layer is again 
used as a mask for the RIE etching of the NbN film. 

III. Experimental and Numerical Results 

All measurements presented in this paper are carried out in 
a fully immersed sample in liquid helium (4.2K). The exper- 
imental setup, used for reflection measurements, is schemat- 
ically depicted in Fig. [2] The samples's RF feedline is con- 
nected to a vector network analyzer using a semi-rigid coax 
cable. The dc feedlines are connected to a dc source-measure 
unit using 4-probe wiring. The laser source has a wavelength 
of 1550nm. 

A. DC I-V Measurements 

The basic I-V characteristics of the HED meander stripline, 
shown in Fig [3] exhibit a highly complex hysteretic behavior. 
Panel [3a) plots nine current measurements, obtained while 
increasing the applied voltage (blue), one on top of the 
other, where each measurement starts at zero voltage and 
ends at a different maximum voltage, slightly above the 
voltages Vcn, depicted in the figure. The corresponding nine 
current measurements, obtained while decreasing the applied 
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Fig. 2. Setup used for reflection measurements. 



voltage are also plotted (red). The magenta and black curves 
correspond to similar measurements, taken while the HED is 
being illuminated. The measurements for low applied voltages 
and currents are enlarged in the insets of Fig |3ja) and (b) 
respectively, where the finite resistance is due to the contacts. 

At panel |3Ja), nine, clearly distinguished, abrupt jumps in 
the measured current are noticed. The number of current jumps 
corresponds to the number of stripline sections that compose 
the meander shape of the HED. Each jump is the result of 
a large increase in the HED resistance due to a transition of 
one section from the superconducting state to the normal one. 
This behavior is typical for a superconducting microbridge 
and is caused by the formation of a hotspot in the bridge area 
[23]. Each critical voltage Vc n i varies, in general, between 
different scans, indicating thus, the stochastic nature of the 
transitions between bistable states. The fluctuation AVcVi in 
Vc-n, between different scans, characterizes the lifetime of 
the pre-jump metastable states of the system. The increase in 
AVcn at high voltages indicates a decrease in the lifetime 
of metastable states because of larger temperature fluctua- 
tions. The combined results of the increasing and decreasing 
applied voltage measurements show that large hysteresis is 
present at all current jumps except for the first one. This 
observation indicates, that only one section at a time can 
be biased into subcritical conditions. Furthermore, only the 
section responsible for the first jump, at V~ci, doesn't suffer 
form hysteresis and thus can repetitively respond to radiation. 
Probably, the cause for this discretization of the critical current 
is the non-uniformity in the meander shape of the HED [23], 
[24]. Our finding clearly shows that the non-uniformity may 
substantially reduces the effective area of the HED, up to a 
fraction of 1/9 of its printed area. 

The same measurements are repeated while constantly il- 
luminating the HED with approximately 27nW IR laser. In 
these measurements the current jumps occur at lower applied 
voltages, V Cn < Vc n - I n addition, AVc n are substantially 
widened. Although the decrease in the critical voltage values 
can be explained by local heating due to the IR illumination, 
the increase in AVc n , especially at low voltages, imply that 
photon absorptions cause discrete events that considerably 
increase the instability of the HED. 

Panel [3b) shows voltage measurements, obtained while in- 
creasing the applied current, with (magenta) and without (blue) 
IR illumination, and while decreasing the applied current, with 
(black) and without (red) IR illumination. Two abrupt voltage 
jumps occur at distinguishable critical currents of Ici — 
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Fig. 3. Basic I-V characteristic of the HED meander stripline. (a) Current 
vs. voltage and (b) voltage vs. current measurements. The magenta (blue) 
curves show an increasing applied voltage and current measurement, with 
(without) laser illumination. The black (red) curves show a decreasing applied 
voltage and current measurement, with (without) laser illumination. The insets 
of (a) and (b) magnify the results for small applied voltages and currents 
respectively. 



4.2/iA and Ic2 — 9.8/iA. IR illumination has a negligible 
measured effect on Ici and a strong effect on Ic2 values. All 
voltage jumps suffer from hysteresis and therefore current bias 
is an unsuitable method for repetitive radiation detection. 

B. DC I-V Effect on the Resonance Lineshape 

Fig. |4] shows several |5n| measurements as a function 
of frequency, in the vicinity of the second resonance mode, 
for various HED resistance values. For clarity, the resonance 
curves are vertically shifted upwards, for increasing resistance 
values. The measurements are obtained while applying vari- 
able voltage Vb, and the resistance is measured simultaneously 
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Fig. 4. Several Sir measurements of the second resonance mode, as a 
function of frequency, for various HED resistance values. The measurements 
are obtained while applying variable voltage and the HED resistance is 
measured using 4-probe technique. Plots are shifted vertically for clarity. The 
inset shows the resistance vs. voltage characteristics of the HED. 
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Fig. 5. Several |Sn| measurements of the second resonance mode, as a 
function of frequency. The measurements are obtained while applying variable 
current. 



with the \Sn\ data using standard 4-probe technique. The RF 
input power is set to — 64.7dBm, where the resonator is in 
the linear regime [25]. The inset of Fig.0]plots the measured 
HED resistance as a function of V&. 

The dependence of the resonance characteristics on the HED 
resistance Rued 7 can be described as followed. At zero applied 
voltage the resonance frequency is fo.Vs — 3.71GHz. At 
very low voltages, as the HED is biased far below critical 
conditions, its resistance is negligible and its influence on 
the resonance curve as well. As the resistance increases, the 
resonance frequency slightly red shifts, and more important, 
the Q-factor i s significantly reduced due to dissipation in the 
HED. This behavior continues up to a point, at i?HED — lfcfi, 
where the resonance curve can be hardly detected. When 
increasing the resistance beyond that point the trend of the 
Q-factor changes, the dissipation decreases, and the resonance 
curve reemerges at a new resonance frequency, fo,v n — 
3.665GHz, red shifted by approximately 45MHz relative to its 
original value. The new resonance Q-factor has a value similar 
to the original one. When further increasing the resistance, the 
trend of the (Mactor continues but no additional resonance 
shift occurs. The behavior of the (Mactor suggests that as 
i? H ED increases, the RF current amplitude of the resonance 
mode in the HED is reduced, due to current redistribution, 
and thus the total power dissipation decreases. 

Similar behavior, with one major exception, can be observed 
under applied current, as shown in Fig|5] The blue, green, red, 
and cyan curves are taken with subcritical 0/iA, 4.33/iA < 
Ici, and over critical 4.39/iA, 5/iA > Ici applied currents, 
respectively, where Ici is the current at which a first jump in 
the measured voltage occurs. There are two well defined res- 
onance frequencies, /o,j s = 3.83GHz, and fo,i n — 3.79GHz, 
which corresponds to applied currents below, and above Ici 
respectively. foj s , and foj n slightly differ from foy s , and 
fo,Vn as the two measurements were taken at different thermal 
cooldown cycles. Low £)-factor curves are absent from this 



measurement because under applied current, the HED can 
not be biased into intermediate resistance values. At bias 
currents below Ici, the HED has low resistance, which only 
slightly increases as the current approaches Ici- As a result, no 
resonance shift occurs, and only the Q-factor slightly reduces 
as the current increases. This behavior changes abruptly once 
the HED resistance crosses a rather low, critical value, Re- A 
self-sustained hotspot is generated [23], quickly expends, and 
the HED becomes resistive. Ici is the bias current at which 
^hed = Rc is obtained. This thermal runaway causes an 
abrupt red shift of A/q = 40MHz in the resonance frequency. 
Further increase of the bias current beyond Ici increases 
the power dissipation and heat generation in the HED. This 
increases the local temperature and dissipation near the HED, 
and thus causes (3-factor reduction. 

C. Resonance Frequency Shift Modeling 

To account for our results we calculate the resonance 
characteristics of our device, as a function of HED resistance. 

As shown in Fig. |6] the ring resonator is modeled as a 
straight transmission line, extending in the ±x directions. The 
HED is represented by a lumped discontinuity, Z = R + jwL, 
connecting x = ±b points together, where R is the resistance, 
uj is the angular frequency, and L is the total inductance of 
the meander shape of the HED. The transmission line has a 
cut at point x = a. The couplings to the RF and dc feedlines 
are neglected. 

The voltage along the resonator's transmission line is given 
by a standing waves expression of the form [26] 

V (x) — { ^ C ° S ~^ ^ Sm —b<x<a 
1 C cos (3 (x — a) + D sin f3 (x — a) a < x < b 

where (3 — 2i:f^fi^/c is the propagation constant along 
the transmission line, / is the frequency, e r is the relative 
dielectric constant, and c is the speed of light in vacuum. 
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Fig. 6. Resonator model. 



The current is given by / (x) = (i/f3Zo) dV/dx, where Zo 
is the characteristic impedance of the line. By applying the 
following boundary conditions: (1) I (a-) — I (a+) = 0, 
(2) I (-b) = I (6), and (3) V (b) — V (-6) = / (b) Z, we 
easily derive a boundary condition equation which yields the 
complex resonance frequencies of the resonator. 

A phase transition of Z from superconducting to normal 
state simultaneously causes changes in its resistive, AR > 0, 
and inductive, AL > 0, parts. Both changes contribute to 
a resonance shift in the same direction. For a very thin 
superconducting films AR ^S> AL resulting in AZ = AR. 
Three fitting parameters are used in the model. Best fit results 
are obtained for R = 0O, Zq = 55S1, and ojL = 1A.5Zq, 
which leads to L = 37.3nH for the second resonance mode. 
R <SC Zq in the superconducting state, and thus negligible. 
The characteristic impedance Zq value is in a very good 
agreement with the designed value of 50S7. The calculated 
inductance L of the meander line at 4.2K without any applied 
current is 6.05nH [27], but this value is strongly dependant on 
temperature and current density, so the fit parameter value is 
in reasonable agreement with the calculated one [27], [28]. 

TablellTlsumm arizes the first three solutions of the boundary 
condition equation for the two extreme cases of AR = and 
AR — > oo. Comparing these results to the measured results 
(taken at two different cooldown cycles), also summarized in 
table [n] shows a good agreement, especially for the second 
and third modes, for which the resonator is designed. 

Fig.Qshows the second resonance characteristics, resonance 
frequency and unloaded damping rate 72 [2], of the experimen- 
tal data (blue) and the numerically calculated data (red). The 
rate 72 is extracted from the data plotted in Fig. fusing the 
method presented in the appendix. The upper subplot shows 
the resonance frequency as a function of AR. The calculated 
resonance frequency, at zero resistance, is /o = 3.791GHz, 
which equals the mean value of the resonance frequency 
measured at different cooldown cycles. In this subplot, the 
calculated data is corrected by — 80MHz to overlap between 
the first calculated and measured point. Both curves show 
the same dependence on AR. The lower subplot shows the 
unloaded damping rate 72, as a function of AR. Also in this 
case, a good agreement with the experiment is obtained, and as 
expected, the measured damping rate exceeds the calculated 
one, due to losses, which are not taken into account in the 
model. 
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Fig. 7. Resonance frequency and unloaded damping rate 72 of the second 
resonance mode. 
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Fig. 8. Normalized voltage amplitudes as a function of the ring's angular 
location, for the (a) second and (b) third resonance modes. 



The coupling between the different modes and the HED can 
be characterized by the current amplitude through the HED. 
The model predicts normalized current amplitudes of 7.3% 
and 5.7% inside the lumped element, for the second and third 
modes respectively. This rather weak coupling is the result of 
the rather high kinetic inductance of the HED. To estimate the 
coupling of modes two and three to the feedline [29], we show 
in Fig. [S] the normalized voltage amplitudes, as a function of 
the ring's angular location. The calculated normalized voltage 
amplitudes, at the feedline coupling location, are 71% and 
92%, respectively. The voltage amplitudes distribution have, 
in general, a strong dependence on the resonance frequency, 
and hence on AR, but because of the rather small resonance 
shift, the voltage amplitudes at the coupling location change 
by less than 2%. 

D. IR Illumination Effect on the Resonance Lineshape 

Fig. [9] plots l^n I measurements with (red) and without 
(blue) IR illumination. The effective IR illumination power, 
impinging on the HED, is approximately 27nW. The RF 
input power is set to — 64.7dBm and the HED is biased with a 
subcritical dc current of 4.14/iA, which only weakly influences 
the resonance curve. When the illumination is turned on, 
the resonance frequency abruptly shifts to a lower frequency. 
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TABLE II 

Resonance Frequency Characteristics 





Numerical Results 


Experimental Results 


n 


fo [GHz] 
AR = 


fo [GHz] 
AR -» oo 


A/ [MHz] 


fo [GHz] 
AR = 


fo [GHz] 
A.R -» oo 


A/ [MHz] 


1 


1.913 


1.873 


39.8 


1.59 


1.58 


10 


2 


3.791 


3.747 


43.7 


3.874, 3.711 


3.829. 3.668 


45, 43 


3 


5.654 


5.62 


34.1 


5.634, 5.38 


5.608, 5.35 


26, 30 
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Fig. 9. Several \Sn\ measurements of the second resonance mode, as 
a function of frequency, while applying sub-critical current, with (red) and 
without (blue) IR illumination. 



The new resonance lineshape has the same characteristics 
as the resonance lineshape measured without illumination 
under supercritical bias current of I = 4.39/iA > Ici- This 
measurement clearly shows that the resonance frequency is 
sensitive to IR illumination. The measured results in this 
experiment yield £Q = 4.14. 

E. Fast Modulation of the Resonance Frequency 

Fast modulation of the resonance frequency is performed 
using the experimental setup depicted in Fig^H The resonator 
is excited by a CW pump signal, at frequency fo = 3.71GHz, 
which coincides with the second resonance frequency. The 
optical signal is modulated at frequency A/, using a Mach- 
Zener modulator driven by a second CW signal, phase locked 
with the first one. The reflected power is mainly composed of 
three tones. One is the reflected pump signal at frequency fo. 
The other two are sidebands, produced by mixing the pump 
signal and the optical modulation signal, and are found at 
frequencies /1.2 = fo ± A/. Occasionally, higher orders of 
the mixed signals are also detected. The amplified reflected 
power is measured using a spectrum analyzer, which tracks the 
fi = fo + A/ tone. No dc bias is needed in this measurement 
as the RF probe signal also serves as a bias signal for the HED. 
This bias scheme has two major advantages over the dc bias 




Fig. 10. Setup used for frequency modulation. 



scheme; first the RF pump bias signal has lower noise, as the 
l/f and line noises are avoided. Second, the RF bias signal 
introduces a strong non-linear mechanism [25], [30], which 
will be discussed in a future publication [31], and produces a 
high internal gain of the induced optical signal. 

The device is not designed for radiation detection. Nev- 
ertheless, we find it useful to characterize the response to 
optical modulation by its noise equivalent power (NEP). 
Fig. ^2 shows the NEP of the device for various optical 
modulation frequencies. Each NEP data point is derived out 
of several reflected power measurements in the vicinity of 
fx, where each measurement is performed with a different 
pump power, ranges between —54.7 dBm and —45.7 dBm. 
The best NEP is measured for a modulation frequency of 
A/ = 10MHz and equals 38fW/\/Hz. A rather low NEP 
of 370fW /vTlz is measured for a modulation frequency of 
A/ = 1.699GHz. For this modulation frequency, f\ coincides 
with the third resonance frequency of the resonator. Higher 
modulation frequencies than the ones presented in Fig. II II are 
also observed, but the measured NEP is relatively poor, mainly 
due to the bandwidth limitations of the microwave components 
composing the experimental setup. 

IV. Discussion 

Fast modulation of the resonance frequency is experimen- 
tally demonstrated. Furthermore, the parametric gain thresh- 
old condition is achieved in a CW measurement. The main 
problem, that currently prevents parametric gain to occur, is 
the relatively low photon flux that impinges the HED. Due to 
losses along the optical path, especially the expansion of the 
Gaussian beam from the tip of the fiber to the HED, the largest 



7 




0.5 1 1.5 2 2.5 3 3.5 
Optical Modulation Frequency [GHz] 

Fig. 11. NEP for various optical modulation frequencies. 



photon flux we currently manage to apply is approximately 
13 photons per modulation cycle, at twice the resonance 
frequency. Taking into account the effective area of the HED 
and its quantum efficiency [32] we estimate that the optical 
power flux is two orders of magnitude lower than the threshold 
power. Therefore parametric gain can not be achieved using 
the present device. Future devices will address this problem. 

V. Conclusion 

We have reviewed the design, fabrication, and measurement 
results of our optically tunable resonator. Resonance frequency 
modulation frequency as high as 4.2GHz is achieved, and a 
NEP of 38fW/\/Hz in the IR is measured. The parametric 
gain threshold condition is achieved in a CW measurement. 
Moreover, the results are shown to be in a good agreement 
with a theoretical modeling. The approach used in this device 
has a great potential of achieving a parametric amplification 
in superconducting resonators. 

Appendix 
Damping Rare Extraction 

The universal expression for the reflection amplitude of a 
linear resonator near resonance is [2], [33] 



S n 



in + (71 - 72) 



(i) 



(71 + 72)' 

where ft = uj p — ujq is the pump angular frequency u> p , relative 
to the angular resonance frequency ujq, 71 is the coupling 
constant between the resonator and the feedline, and 72 is the 
unloaded damping rate of the resonance. The damping rate is 
numerically extracted by expanding Eq. Q to first order in ft, 

Sn =r + r 1 n + O(n 2 ) , 
where ro is the Sn value at the resonance frequency, 

72 - 7i 



7i + 72 ' 



and r\ is the slope of the imaginary part of Sn, 

n = — 1 2- 

(71 +72) 

Note that the extraction of r\ is less accurate for low £>-factor 
curves, and thus the calculated loss factor suffers from a rather 
large impreciseness at that regime. 
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